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ABSTRACT 

The bulk composition of an exoplanet is commonly inferred from its average density. For small 
planets, however, the average density is not unique within the range of compositions. Variations 
of a number of important planetary parameters — which are difficult or impossible to constrain from 
measurements alone — produce planets with the same average densities but widely varying bulk com- 
positions. We find that adding a gas envelope equivalent to 0.1%-10% of the mass of a solid planet 
causes the radius to increase 5-60% above its gas-free value. A planet with a given mass and radius 
might have substantial water ice content (a so-called ocean planet) or alternatively a large rocky-iron 
core and some H and/or He. For example, a wide variety of compositions can explain the observed 
radius of GJ 436b, although all models require some H/He. We conclude that the identification of 
water worlds based on the mass-radius relationship alone is impossible unless a significant gas layer 
can be ruled out by other means. 
Subject headings: extrasolar planets 



1. INTRODUCTION 

Out of over 250 exoplanets known to date, over 20 
are known to transit their stars. Transiting planets are 
important because we can derive the precise mass and 
radius, and can begin to determine other planetary prop- 
erties, such as the bulk composition. 

Much attention has been given to "ocean planets" or 
"water worlds" , planets composed mostly of solid water 
(jKuchnerl [20031 : iLeger et al J (2004( 1 . If a water world is 
found close to a star, it would be strong evidence for mi- 
gration because insufficient volatiles exist near the star 
for in situ formation. The proposed identification of wa- 
ter worlds is through transits. From a measured mass 
and radius a low-density water planet could potentially 
be identified. 

We examine the possibility that water worlds can- 
not be uniquely identified based on the mass and ra- 
dius of a transiting planet. An alternative interpretation 
could be a rocky planet with a thick hydrogen-rich at- 
mosphere. Most authors have assumed that solid plan- 
ets in the 5 t o 10 M ffl range have a n insignificant amount 
of hydrogen ([Valencia et al.|[2006allb]: iFortnev et al.ll2007l : 
ISeager et aTll2007t ISelsis et al.ll2007t ISotin et alJl2007t l. 

Exoplanets have, however, contradicted our basic as- 
sumptions before. Notable examples include: the exis- 
tence of hot Jupiters; the predominance of giant plan- 
ets in eccentric orbits; and the gas-rock hybrid planet 
HD 14902 6b with its ~ 60 M ffl core and ~ 30M e H/He 
envelope ((Sato et al.ll2005h . 

We adopt the idea that a wide range of atmospheric 
formation and loss mechanisms exist and can lead to 
a range of atmosphere masses on different exoplanets. 
We explore the mass-radius relationship for the lowest- 
mass exoplanets yet detected (~ 5-20 M®) in order 
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to identify potential ambiguities that result from the 
presence of a massive atmosphere. We explore atmo- 
spheres ranging from ~ 10~ 3 M ffi (10 x Venus' atmo- 
spheric mass) to ~ l Mffi (the estimated mass of Uranus ' 
and Neptune's H/He (|GuillotJl2005l: iPodolak et alJfl995t 
iHubbard et al.l(l99l( l), with a focus on the smaller mass 
range. We also explore potential composi tions for the 
transiting Neptune-size planet GJ 436b ((Butler et al.l 
l2lMlGlllon et al.ll2007bh . 

2. MODELS 

Our model assumes a spherical planet in hydrostatic 
equilibrium, with concentric shells of different compo- 
sition sorted by density. We solve for the mass m(r), 
pressure P(r), density p(r), and temperature T(r) using 
the equation for mass of a spherical shell, 



dm(r) 
dr 



Airr 2 p(r), 



the equation of hydrostatic equilibrium, 

dP(r) _ Gm(r)p(r) 
dr r 2 
and the equation of state of the material. 



(1) 



(2) 



p(r) = f(T(r),P(r)). 



(3) 

To build a model planet, the composition and mass 
of up to three layers are specified. Starting from a cen- 
tral pressure adjusted to achieve the desired mass, the 
differential equations for P(r) and m(r) are numerically 
integrated outward, and p(r) is calculated from the equa- 
tion of state for the material. At the boundary between 
layers, P(r) and m(r) are used as the initial conditions 
for the next layer. The integration is stopped at the 
outer boundary condition P(r) — 1 bar, approximately 
where the atmosphere ceases to be opaque at visible and 
infrared wavelengths 3 . 

3 We ignore any wavelength-dependent effects on potential mea- 
sured radii of extrasolar planets. 
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The temperature of a planet's gas layer has a sig- 
nificant effect on the radius. Instead of computing a 
cooling history to obtain T(r), we calculate T(r), sep- 
arately in three different regimes. Within the solid por- 
tion of the planet (r > r so ud), temperature has little 
effect on the final radius, an d is assumed to be isother- 
mal (see i Seager et al.l 120071 ). In the deep hydrogen- 
helium layer where the pressure is greater than 1 kbar 
{r solid < t < rikbar), we assume that convection domi- 
nates and the temperature follows an adiabat tied to the 
entropy, S, at 1 kbar. At lower pressures (r > rikbar), 
we use the radiat ive eq uilib rium gra y analy tical model of 
IChevallier et al.l (|2007f l and lHansenl (|2007f ) for irradiated 
atmospheres, 



1981. and H 2 Q are described in more detail in 



Seager et all (|2007f ). 



T(r) 



T eq F ( T ' MO, 7), 



(4) 



where r is the optical depth, /iq = cos 9q where 9q is the 
angle of incident radiation with respect to the surface 
normal, and 7 is a parameter that accounts for the al- 
titude at which radiation is absorbed. We convert from 
optical depth to pressure through a constant scaling re- 
lation under the assumption of hydrosta tic equil i brium 
and constant opacity per gram, following Ha nsen! (|2007l ) 
but taking into account the planet's surface gravity. Here 
T e ff is the effective temperature in the absence of stel- 
lar irradiation, representing energy from internal sources. 
T eq is the equilibrium temperature which represents heat- 
ing from the parent star. 

We emphasize that by choosing an effective 
temperature — or adiabat — we are subsuming a cooling 
calculation. This is common practice for modeling the 
interior structu re of solar sys tem giant planets (e.g. , 
Stevensodll982tlHubbard et al.lfl99lt iMarlev et~aDl995t 



Saumon. fe Guillotl 12004) and was also the case even be 



fore their T e ffS and gravitational moments were known 
or w ell known (jDemarcusI 119581 : IZapolskv fe Salpeterl 
1969). In fact, even with a known T e tf cooling 
calculations only match observed parameters for 
!.g.lHubbard|[l977l): Satur n is hotter than ex- 



pected dPollack et al.l 119771 : iStevenson fc Salpeterl 119771 : 



Jupiter (e.p 
(IPol 

iFortnev fc Hubbard! l2003f) and Uranus a n d Neptune are 
colde r than expecte d dStevensonl 11982b iPodolak et al.l 
119911 : iGuillotl l2005h . An additional motivation to 
adopt a simple framework for temperature is the huge 
uncertainty in interior parameters for rocky planets. 
This includes the temperature-dependent equation of 
state of liquids and solids at high temperature and 
pressure, the temperature-dependent viscosity, and 
the effect of tides on the cooling history for eccentric 
exoplanets. We further note that our choice of radiative 
equilibrium down to 1 kbar is b ased on the ir r adiate d 
atmosphere/interior models by IFortnev et al.l (|2007f ). 
For our fiducial model we choose T e ff to be similar to 
Earth's and Uranus's. While the largest uncertainty 
in our treatment is the qualitative choice of T e ff, we 
subsequently vary over a reasonable range of effective 
temperatures (or adiabats). 



We use the H/He equation of state from lSaumon et al.1 
(fl995h . ignoring the "plasma-phase transition" which 
may be a numerical artifact (D. Saumon 2007, private 
communicati on). The equations o f state for the solid 
materials Fe ([Anderson et al.ll2001h . MgSiOs perovskite 
(iKarki et all |2000|), (Mg^Fe)Si0 3 dKnittle fc Jeanlod 



3. RESULTS AND DISCUSSION 

Our fiducial planet consists of a 30% Fe core and a 70% 
MgSiC>3 mantle, roughly analogous to Earth. We used 
a H/Hc mixture with helium mass fraction Y = 0.28 
(the He mass fraction of the solar nebula). We chose 
T eq = 300 K, based on the observation that a planet 
around an M-dwarf at an orbital distance of 0.1 AU has 
a similar equilibrium temperature to Earth (assuming 
similar albedos). We set T e ff — 30 K, similar to Earth 
and Uranus . For the atmospheric parameters, we fixed 
fio = cos 60° and 7 = 0.1 to represent radiation absorbed 
deep in the atmosphere 5 . We later investigate variations 
on Y", T e ff, T eq and 7. Figure 1 shows a plot of the 
mass-radius relationship for fiducial planets of masses 5, 
10, 15, and 20 Mq. For each planet mass, we added 
atmospheres ranging in mass from 0.001-1 M@. 

A robust finding for all models is that a small amount 
of gas creates a large radius increase. While this result is 
expected, the radius increase is far more dramatic than 
anticipated. 

For example, an H/He atmosphere of ~ 0.001 by 
mass — only ten times greater than Venus' atmospheric 
mass fraction — is required for a noticeable radius in- 
crease. As seen in Figure 1, adding a hydrogen-helium 
atmosphere with just 0.1% of the mass of a 10 M© rocky 
planet results in a 5% increase in the planetary radius — 
within a measurement precision that has been obtained 
for currently known transiting planets. 

As a second example, adding a gas layer of H/He equal 
to 1% of the mass of our fiducial planets increases the 
radius by ~ 20% of the original planet radius, or by about 
0.35 i?® for the planet masses we considered. 

Our major finding is that that exoplanets with a sig- 
nificant H/He layer cannot be distinguished from water 
worlds, based on M p and R p alone. For our fiducial solid 
exoplanets, adding up to 5% H/He by mass (for 10 M e 
planets) is sufficient to push the planet's radius through 
the entire range of radii corresponding to solid planets 
with no gas, including planets with up to 100 percent 
water composition. While we have not completed an ex- 
haustive study of possible compositions, we find the non- 
uniqueness of water planets to be valid for any conditions 
we investigated. 

This generic finding holds for a wide range of as- 
sumptions of assumed temperatures. Taking our fiducial 
model, we vary T eq , T e ff, and 7 individually. For a 10 
M© solid planet with an additional 0.1 M© H/He atmo- 
sphere, increasing T eq from 300 to 500 K increases the 
radius by about 1% (Figure 2). For the same planet 
varying T e ff from 10 K to 50 K results in an 8% increase 
in radius. While large, this value is comp arable to ex- 
pected r adius uncertainties fo r these planets (|Gillon et al.1 
l2007blial : lDeming et al.ll2007h . Varying the altitude where 
radiation is absorbed (specified by 7) has a much smaller 
effect on the planet radius. Varying 7 from 0.1 to 10 
causes the radius to decrease by 0.2%. 



4 Earth has 44 X 10 12 W (Pollack ct al. 1993) and Uranus has 
340 X 10 12 W of energy flow llPearl et al.|[T990T) . 

5 In comparison 7 = 10 would correspond to absorption high in 
the atmosphere. 
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Fig. 1. — The increase in radius due to adding H/Hc to a solid 
planet. A H/He layer of 0.002-1 M® is added to a solid planet of 5, 
10, 15, or 20 M®, with fiducial model parameters (30% Fe and 70% 
MgSi0 3 ). The black points are for atmospheres at 0.01 Af® and 
every 0.1 M® afterwards. The mass-radius relationship of solid 
planets with no gas is plotted for comparison. The water (blue), 
rock (red) , and iron (green) curves are taken from ISeager et al.l 
(2007) and represent homogeneous solid planets. Intermediate 
compositions for differentiated planets are, from top down: dashed- 
blue, 75% H 2 0, 22% MgSi0 3 , and 3% Fe; dashed-dotted-blue, 48% 
H 2 0, 48.5% MgSi0 3 , and 6.5% Fe; dotted-blue, 25% H 2 0, 52.5% 
MgSi0 3 , and 22.5% Fe; dashed-red, 67.5% MgSi0 3 and 32.5% Fe; 
and dotted-red, 30% MgSi0 3 and 70% Fe. In general, the addition 
of a gas layer of up to ~ 5% of the solid planet mass will inflate 
the radius of a rocky-iron planet through the range of radii cor- 
responding to water planets with different water mass fractions. 



Teq (300. 400 and 500 K) 




10. 10.1 10.2 10.3 10.4 10.5 



TeffdO, 30 and 50 K) 




Fig. 2. — The effects on the radius of varying the equilibrium 
temperature (top) and effective temperature (bottom) for a 10 Af® 
planet with otherwise fiducial parameters. T eq values of 300, 400, 
and 500 K are plotted to simulate the effect of uncertainty in the 
orbital parameters and albedo on the expected radius. T e f f values 
of 10, 30, and 50 K are plotted on the same scale as the T eq plot, 
to show uncertainties in the planet's interior temperature. Uncer- 
tainties in the internal energy of a planet lead to large variations 
in radii for a given mass, showing how temperature is a large un- 
certainty in the interpretation of a planet's internal composition. 



A corollary of our main result is that when a planet has 
a significant H/He atmosphere there is a wide degener- 
acy in allowable internal composition. This is not just 
compositional, but also relates to the trade off of tem- 
perature and mass of H/He gas. ft could be argued that 
specifying a planet's composition implies a particular in- 
ternal thermal profile derived from a consistent cooling 
history. As addressed in §2, the many unknowns and 
free input parameters for rocky planet interiors — such as 
the possible differences between atmospheric and interior 
compositions, equation of states, and the effect of tides 
on the planet's cooling history — prevent a self-consistent 
solution for the present time. 

How could a 5-20 M ffi exoplanet get a substantial 
H/He layer? Two different scenarios may produce them: 
direct capture of gas from the protoplanetary disk (possi- 
bly modified by the escape of some fraction of the original 
gas) or outgassing during accretion. 

A planet may capture and retain up to f to 2 M® of 
H/He if the planetary core did not grow quickly enough 
to capture more before the gas in the disk evaporated 
(as is the paradigm for Uranus and Neptune). Alterna- 
tively, for short-period exoplanets, a f to 2 Af® H/He 
envelope may result after substantial loss of an initially 
massive gas envelo p e from irradiated e vaporation (e.g. 
iBaraffe et al.ll2006T ) . lAlibert et alJ <|2006l ) consider atmo- 
spheric evaporation during migration, and conclude that 
the f0M e innermost planet in HD 69830, at 0.08 AU, 
kept ^2 of H/He over the 4 Gyr lifetime of the star. 

Little attention has been given to the mass and compo- 
sition of exoplanet atmospheres from outgassing. Venus' 
atmosphere is 10 _4 M e ; if Venus had a surface gravity 
high enough to prevent H escape, its atmosphere would 
be over 10 _3 M®. Even more massive H-rich atmospheres 
are possible. If a massive iron-silicate planet formed with 
enough water, the iron may react with the water dur- 
ing differentiation, liberating hydrogen gas (jRingwoodl 
Il979t IWaenke fe Dreibusl 11994ft . L. Elkins-Tanton et al. 
(in prep.) estimate that the maximum H component is 
about six percent by mass for a terrestrial-composition 
planet. For a fOAf® planet this would result in a 0.6A1® 
H envelope. 

For short-period, low-mass planets, theoretical argu- 
ments of atmospheric escape may be the best way to 
identify a water w orld based on the mass and radius mea - 
surements alone (jLeger et al.l |2004 ISelsis et all [2007). 
Indeed, our assumption of H/Hc atmospheres for exo- 
planets relies on the condition that atmospheric mass loss 
has not evaporated all of the H/He. In the absence of 
hydrodynamic escape, the exospheric temperature (and 
not the atmospheric T e ft) drives the thermal Jeans es- 
cape of light gases. Eart h and Jupiter both have ex obase 
temperatures of 1000 K (|de Pater fe Lissauerj|200lD . sig- 
nifican tly above t heir T e /f of 255 K and 124.4 K respec- 
tively (|Coxll2000| ). Uranus and Neptune have exobas e 
temperatures around 750 K (|de Pater fe Lissaueij|2001h . 
We note that because GJ 436 must have at least 1 M® 
of H/He, its exospheric temperature is not too high. On 
the other extreme, planets of 5M® would require very low 
exospheric temperatures (~ 300 K) to retain a massive 
atmosphere over the course of billions of years. Never- 
theless, a young 5Af® Earth-mass planet with a captured 
atmosphere could still have a H/He atmosphere and an 
old 5M® planet could retain a substantial He fraction, 
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Fig. 3. — Density vs. radius for three different potential compo- 
sitions of GJ436b. From top to bottom: solid (black) curve, 19.0 
M e core (30% Fe, 70% MgSi0 3 ) with 3.2 M ffi H/He (Y=0.28); 
dotted (red) curve, 20.0 M® core (100% (Mg,Fe)Si0 3 ) with 2.2 
M e H (Y=0); dashed (blue) curve, 20.5 M e core (90% H 2 0, 10% 
MgSi0 3 ) with 1.7 M e H/He (Y=0.28). All three planets have the 
same total radius (4.3 Rqj) and total mass (22.2 Mj). 

making its compositional identification ambiguous. 

We now turn to a qualitative study of GJ436b, to show 
that the interpretation that GJ436b is a water world 
akin to Uranus and Neptune (jGillon et al.ll2007a|) is not 
the only possibility. We consider the G J 436b values 
M P = 22.2M© and R p = 4.3i?© from iDeming etldl 
(2007). The internal structure in Figure 3 shows how 
three planets with very different internal compositions 
can have the same total mass and radius. We first ex- 
plore a planet similar to our fiducial model: a 22.2 M® 
solid planet with Earth-like iron/rock mass ratio (30/70), 
T e ff = 30 K, and T eq — 600 K in rough agreement 
with the orbital parameters (assuming an albedo of 0.1). 
By adding ~3.2Af® of hydrogen-helium to the 19.0 Af® 
solid planet, we are able to reproduce GJ 436b's radius. 
We note that the mass of gas is 15% of the solid mass, 
likely too much to have originated from outgassing, and 



so capture must be at least partially invoked to explain 
such a massive atmosphere. The second composition for 
GJ 436b we considered is for water worlds, one with 
a 50% water mantle (by mass) and 50% silicate core, 
and another with 90% water mantle and a 10% silicate 
core. These planets also need some H/He to match the 
known radius, 12% and 8% by mass, respectively. The 
third model approximates planets with atmospheres cre- 
ated from outgassing, considering an extreme scenario 
where all of the available water has oxidized iron, leav- 
ing a 100% (Mg,Fe)Si03 solid planet core. To match 
the observed radius, a 22.2 M® planet requires ~ 2.2M® 
of H alone, a case that assumes no initial trapping and 
subsequent outgassing of He. The pure-hydrogen atmo- 
sphere is 10% of the mass of the solid planet. This is 
above the theoretical maximum of outgassing based on 
observed abundances of metallic iron in chondritic mete- 
orites from our solar system (see Elkins-Tanton et al. in 
prep). Although not an exhaustive study, the range of in- 
terior compositions illustrates the variety of possibilities, 
though all models require some H/He. 

While our study is preliminary, we make the robust 
point that H-rich thick atmospheres will confuse the in- 
terpretation of planets based on a measured mass and 
radius. This point is independent of the uncertainties re- 
tained by our model including T eq , T e ff, the mass frac- 
tion of H/Hc, and the mixing ratio of H and He. We 
find that the identification of water worlds based on the 
mass-radius relationship alone is impossible unless a sig- 
nificant gas layer can be ruled out by other means. Spec- 
troscopy is the most likely means and may become rou- 
tine with transit transmission and emission spectroscopy, 
and eventually with spectroscopy by direct imaging. 



We thank Mark Marley, Jonathan Fortney, and Wade 
Henning for useful discussions. 
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